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We have postulated a single pole for the self-energy, E(k,u), looking for the consequences on the 
one-particle Green function, G(k, uj) in the Hubbard model. We find that G(k, uj) satisfies the first 
two sum rules or moments of Nolting (Z. Physik 225, 25 (1972)) for any values of the two unknown 
k parameters of E(k, uj). In order to find these two parameters we have used the third and four sum 
rules of Nolting. G(k,ui) turns out to be identical to the one of Nolting (Z. Physik 225, 25 (1972)), 
which is beyond a Hubbard-I solution. Furthermore, we have also postulated a two-pole Ansatz for 
E(k,oj) which is equivalent to a three-pole Ansatz for G(k,u). In the latter case, we have fixed 
two of the frequencies. This Green's function is also beyond the Hubbard-I solution. We present 
numerical results for the band structure and the spectral weights. 
Pacs numbers: 74.20.-Fg, 74.I0.-z, 74.60.-w, 74.72.-h 



After the discovery of the high-T c materials jj], the study of correlations has gained interested due to the fact 
that there is the belief Q that the normal properties of these materials could be explained in the framework of the 
Hubbard model Q, since electron correlations are strong, i.e., the on-site electron-electron repulsions U are much 
larger than the energies associated with the hybridization of atomic orbitals belonging to different atoms This 
Hamiltonian is a kind of minimum model [^) which takes into account quantum mechanical motion of electrons in a 
solid, and nonlinear repulsion between electrons. Even though this model is too simple to describe solids faithfully, 
serious theoretical studies have revealed that to understand its various properties is a very difficult task, since is the 
simplest many-body Hamiltonian one can write down and which cannot be reduced to a single-particle theory Q . Its 
study will prove useful in developing various notions and techniques in statistical physics of many degrees of freedom. 
Besides cuprate superconductors, the Hubbard model may also be applicable to describing the metal - insulating 
transition in materials like La x Sri- x TiO^ and V^-yOs in which paramagnetic metal, antiferromagnetic insulator, 
and a phase antiferromagnetic metal are clearly separated in temperature - pressure phase diagram §. All these 
features, found in the Hubbard model, make it suitable to reproduce experimental data. Additionally, the study of 
correlations by itself in the Hubbard model is a rewarding task since it will shed light on still unsolved points of the 
novel materials. For example, at high temperatures (T c 30 — 130 K) these HTSC cuprates, which are poor conductors, 
become superconductors. This feature is strange indeed because the Coulomb repulsion is strong. Furthermore, the 
behavior of these materials at T > T c is even more puzzling than the superconductivity itself. Contrary to the 
predictions of the Fermi liquid theory, the resistivity at T > T c and optimum doping is linear in temperature, i.e., 
R w T . This suggests a very strong scattering of elementary excitations. A discussion of the possible breakdown of 
Fermi liquid theory is given in Ref . . 

In this paper, we will adopt a given pole structure in the self-energy and see its effect on the one-particle Green 
function. This is done with the idea of reducing the amount of computational effort which is implicit in the moment 
approach of Nolting which is normally constrained to two poles in G(k,w). At the same time, without additional 
effort, we monitor the self-energy which is an important quantity for life-time considerations. We find that by using 
this approach (one pole in S(k,o;) and the sum rules for the spectral functions), the one-particle Green function 
goes beyond the Hubbard-I solution, i.e., it identically satisfies the first two sum rules and its parameters have to be 
calculated self-consistently. This has motivated us to consider two poles in the self-energy and to see the effects on 
G(k, The general result that we gain is that iV-poles in £(k, uj) is equivalent to N + 1-poles in G(k, uj). Along the 
way, we discuss the possibility of obtaining a third peak in G(k, uj) by means of this technique. 

The model we study is the Hubbard model [|| 

H = t^^cl c-r H n^- — ucl o , (1) 

where ct (c* ) are creation (annihilation) electron operators with spin cr. = cl C| . U is the local interaction, 
\x the chemical potential and we work in the grand canonical ensemble. We have adopted Einstein convention for 
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repeated indices, i.e., for the N s sites i, the z nearest-neighbor sites and for spin up and down [a = ±1). tfj = —t, 
for n.n. and zero otherwise. 

The one-particle Green function, G(k, ui), is expressed in terms of the self-energy, E(k, w), as 



G(k;w) = 



1 



- £k - E(k, 



(2) 



where e(k) = — 2t(cos(k x a) + cos(k y a)), and £k = e(k) — /j,. 
We adopt the following Ansatz for S(k, u>): 

E(k,w) = + 



g(k) 

j — Ol 



(3) 



As S(k, a;) has dimensions of energy, the still unknown parameter a^ has dimensions of (energy) 2 . a(k) is ] kind 
of a spectral weight and Ok is the energy spectrum of the self-energy. We will calculate a(k) and ilk- Along the way, 
we will give a simple physical interpretation for a(k). We have included the Hartree shift directly in the self-energy. 
In Eq. (□), p is the carrier concentration per spin orientation, i.e., p = n/2. We are assuming that we are in the 
paramagnetic phase 

By using Eq. (ph into Eq. (0), we get that the one-particle Green function has two poles. It can be written as 



G(k,w) 



«i(k) 



« 2 (k) 



to — wi(k) w — w 2 (k) 
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w 2 (k) - wi(k) 

From Eqs. (^) we inmediately see that the following sum rules or moments |l(J are identically satisfied: 

ai(k) + a 2 (k) = 1 , 
ax (k)wi (k) + a 2 (k) u> 2 (k) = £ k , 
Ck = e(k)+pf/ . 



(5) 



(6) 



Eqs. (0) are the first two sum rules for the spectral functions of Nolting fTcfl . In order to evaluate a(k) and Ok, we 
use the next two sum rules of Nolting. This gives: 



^ 2 (k) ai (k) + w 2 2 (k)a 2 (k) = £ + a(k) = oa(k) , 
w?(k) ai (k) + w 3 (k)a 2 (k) = C k + (fik + 2? k )«(k) = o a (k) , 

where a 2 (k), a 3 (k) are given in Ref. M]- Q as 

Oa(k) = e 2 (k) + 2p£/e(k) + pU 2 
03 (k) = e 3 (k) + 3t/e 2 (k)(2 + p)pC/ 2 e(k) + p(l - p)ELB(k) + pt/ 3 



where -B(k), in the spherical approximation of Nolting, is given by 
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p(i- P )B(k) = Y,J2 



3=1 k 



( Wi (k) - e(k)) - 1 



; (k)/K(k)) £ (k) 



(7) 



(8) 



(9) 



where /(x) is the Fermi distribution function. By solving Eqs. ([?]) we find: 
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a(k) = P (l - P )U 2 , 

n k = {i- P )u + B{k) , (io) 

B(k) is the narrowing band parameter denned by Nolting. This narrowing band parameter has to be calculated 
self-consistently and in the spherical approximation of Nolting is k independent. By combining Eqs. (0J^) we find 



1 



2 

where 



Wl (k) = ^ H(k) + [K(k)}^ ; w 2 (k) = ~ ff(k) - [/if(k)]V=« , (11) 
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HQs) = e(k)J7 + B(k) , 
J fC(k)EE(e(k)-C/-B(k)) 2 + 4 i o[/(e(k)-B(k)) . (12) 

Eqs. ( |l2|) are nothing that the solutions given by Nolting Jic[| by solving the four unknown ttj(k) and Wj(k), with 
i = 1,2. We have regained Nolting's solutions in a much easier way, starting from the self-energy, while Nolting does it 
from the Green function itself. At the same time we have been able to identify the self-energy for the two-pole ansatz 
for the one-particle Green function. Even more, we see that the self-energy has the meaning of being an expansion in 
powers of U . In our case, due to the choice of a single pole for the self-energy, a(k) is of second order in U. Remember 
that the first order expansion has been taken into account by the Hartree shift. Then, if we took two poles for the 
self-energy, then a(k) has to go up to order U 3 . This is a conjeture we will justify later. We would like to point out 
that in the atomic limit, i.e., t = 0, we recover the two Hubbard bands, i.e., 

Wi(k) = -/i + {7 ; w 2 (k) = - p , 

«i(k)=p ; a 2 (k) = l-p , (13) 

as it should be. As our solutions given by Eqs. (|l^) satisfy the first four sum rules or moments, our Green function 
is beyond the Hubbard-I solution [Q and as a consequence we have an improved solution pT| . The drawback of 
the Hubbard-I solution, i.e., a gap for any value of the interaction, was pointed out by Laura Roth |r| many years 
ago. This gap exists, for all values of U, in the spherical approximation of Nolting |10[ |- [ pd| . It is known that the 
Hubbard-I solution does not show ferromagnetism. This difficulty is solved by the presence of the jB-term in Nolting 
solution which is spin-dependent. 

Eqs. (ra) contain interesting solutions. One of them we consider next. For example, if we take 



we inmediately get 



and 



« x (k) = - c^ 2 (k) , (14) 

wi(k) = [e k +p(i- P )u 2 ] 1/2 =E k , as) 



a 2 (k)=l- ai (k) = i(l--^) . (16) 

Eqs. dT^| l6|) are the well-known spin-density- wave (SDW) solutions at half-filling, with a gap given by U/2 — 
{—U /2) = U Hl4|] . Thus, our solutions reduce to known results for some values of parameter space. Chubukov and Morr 
Jl5[ have discussed the mean-field SDW result within the effective spin-fermion model. According to these authors a 
one-pole self-energy close to half-filling includes only exchange of longitudinal fluctuations. Substituting this diagram 
(self-energy) into Dyson equation produces inmediately the mean-field SDW result. Chubukov and Morr's analysis 
|jl5f is based on the theory of Kamp and Schrieffer jl6) where it is assumed that the longitudinal spin susceptibility, 
X s (q, u>) has a ^-functional peak at zero frequency and momentum transfer Q = (n, n), i.e., Xs(qj w ) = ^(Q - Q)^ 1 - )- 
Chubukov and Morr's paper studies the Fermi surface evolution with interaction. They find that the SDW mean-field 
solution (valid for large U) yields a small Fermi surface centered at (tt/2,tt/2). Notice, however that they include a 
term t' in the free-energy band. Let's mention that the SDW state has not been observed in FLEX calculations |17| ] 
where short-wavelength spin fluctuations are taken into account. By using the Green function given by Eqs. (pl[L2fr, 
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we can evaluate the dynamical conductance, cr(w), vs lu, away from half-filling |T^ ]. We mention that the two-pole 
one-particle Green function (Eq. (|^) has been used to calcultate the static spin susceptibility, \s{T), pOj]- [jll| and 
the results compare rather well with Quantum Monte Carlo simulations [18| . Another quantity which can be easily 
calculated is the specific heat, cy, vs temperature, T J O | . See also Ref. fll9fl . We would like to comment on the fact 
that our two energy branches (second line of Eq. (|l^)) produce a gap with a value of U . This gap is evidently no 
perturbative contrary to the view of Ref. [[l6j . 

Now, let us propose a two-pole Ansatz for the self-energy, given by 

Z(k,u) = pU+- Jffi . (17) 

The assumption of Eq. (|l^) is equivalent to have a three-pole structure for the one-particle Green function, i.e., 

G(k,.),^^ + ^^ + ^- . (18) 

LU — CiJi(k) CO— W 2 (k) LO — L03{K) 

The Wj(k)'s, j=l,2,3, are the roots of the following cubic equation 



U! 
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AtMcv 2 + A 2 (k)uj - A 3 (k) = , (19) 



where 



A 1 (k) = n 1 (k) + n 2 (k) + & ; A 2 (k) = &(n 1 (k) + n 2 (k)) + f2 1 (k)n 2 (k) , 

A 3 (k) = a(k)+0 1 (k)0 2 (kX k . (20) 



The Oj(k)'s, j=l,2,3, are given by 

«i(k) 
a 2 (k) 
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-n 2 (k)) 
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- Wl (k))(w2(k) 


-wa(k)) 


(w 3 (k) 


- fii(k))(w 3 (k) 


-n 2 (k)) 


Mk) 


-wi(k))(w 3 (k) 


-wa(k)) 



= v^aw "iwjww " 2 y*y; 



Now, as we are after a simple solution to our Eqs. (|19|-pi[), we are going to impose the strong condition on the 
cubic equation (Eq. (|l9|)) that 

^ 3 (k) = . (22) 



Then, Eq. ( |19| ) becomes 

uj 2 - A^lv + A 2 (k) = , 
a(k) = - ni(k)n 2 (k)& . (23) 

If we require that the other two roots are antisymmetric, we get 

j4i(k) = fii(k) + n 2 (k)+& = o , 
^ = -[fii(k)n 2 (k)) + &(fii(k)+n 2 (k))] . (24) 

We would like to point out that this condition (the first line of Eq. fl24|)) is not really necessary. It has been used 
here to keep the algebra as simple as possible. Otherwise, higher order moments will be needed. 
It is easy to express the a.j (k) 's as 



ar(k) = ^2 (O k -Oi(k))(O k -0 2 (k)) , 

a 2 (k) = ^(f\ + fMk))(tt k + tt 2 (k)) , 
»i(k)Q 2 (k) 

«3(k) = ^ • ( 25 ) 
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The first two sum rules or moments are satisfied identically since 



ai(k) +a 2 (k) + a 3 (k) = 1 , 
fi k (a 1 (k)-a 2 (k)) = -(n 1 (k)+n 2 (k))=a . (26) 

The second moment is satisfied due to the fact that Ai(k) = 0. (See Eq. (|23)). Let us point out that the first sum 
rule is always satisfied, without the requirement that ^(k) = 0. 
The third and four moments are given by 

ni ( ai (k) + a 2 (k)) = ill + n!(k)0 2 (k) = a 2 (k) , 
Q|(ax(k) - « 2 (k)) =- fi^(Oi(k)+0 2 (k)) =a 3 (k) . (27) 



Combining Eqs. (|(i]|27|) we end with 



<* = ?i • (28) 
ai(k) 



The other unknowns are easily found. They are: 

-£ k ±K + 4(^-a 2 (k))] 1/2 



n ll2 (k) 



/2 



ai :2 (k) = 



Ck(a 2 (k)±[aa 3 (k)] 1 

2a 3 (k) ' 

x , Cka 2 (k) 

a 3 k = 1 — — , 

«s(k) 

a(k) = a 3 (k)-£ k a 2 (k) . (29) 

Now we are in a position to justify the fact that X(k, uj) represents an expansion in powers of U . (See the discussion 
after Eq. ([l2|)). From the four line of Eq. (^9|) we can calculate the explicit form of the self-energy weight. We use 
the expressions for a 2 (k) and a 3 (k). The result is 

a(k) =p(l- P )U 2 [U + 2e(k) + B(k)} . (30) 

Then, we see from Eq. (|30| ) that we have an expansion in terms of U. Just imagine the atomic limit!. Let us say that 
the solutions we have found (Eqs. @ H^USO 

can also be obtained from the three pole Ansatz for the one-particle 



Green function, if we require that the first four moments are satisfied together with the conditions of a peak at zero 
frequency and the other two peaks at antisymmetric frequencies. In Eq. (pfl) B(k) is given by 



3 r 2 



p(i-#)=EE 

j=l k 



{J v^j(k) -e k ) 



a,-(k)/( Wi (k)) e (k) . (31) 



In Fig. 1 we present the energy spectra, i.e., u>i(k), i = 1,2,3 vs k along the diagonal of the Brillouin zone, i.e., 
k x = k y = k for T/t = 0.5, p = 0.5 and p = U/2, with U/t — 7.0. In Fig. 2 we show the spectral weights, aj(k), 
i = 1,2,3, for k along the diagonal of the Brillouin zone for the same parameters as in Fig. 1. In Fig. 3 we plot 
f2j(k), i = 1,2 along the diagonal of the Brillouin zone for the same parameters as in the first two figures. We observe 
from all these figures that there are regions of k which are forbidden, i.e., for which wi(k become imaginary. We have 
neglected these points. However, they could be taken into account if we include lifetime effects. Our choice of a delta 
function for the spectral function, A(k, lj), forbids us to consider imaginary frequencies. 

In summary we have proposed a self-energy of N poles which implies that the Green function is composed of 
N + 1 poles. This has been accomplished by the use of Dyson's equation (E q. H ). For the case of a one-pole for the 
self-energy we have reproduced with minor effort the Nolting's solution |L0]]- Jllf ! The one-pole self-energy leading to 
Nolting's solution puts on firm grounds the results found in Refs. p"l[ where it was argued that the two branches of 
the one-particle Green function were due to a single branch in the self-energy. For the case of two-pole Ansatz for the 
self-energy we have worked out the simple case of what we call the Kondo peak. In the later case, we have been able 
to solve all the unknown parameters. In both cases, the first two sum rules are satisfied identically. In both cases, 
our Green's function goes beyond the Hubbard-I solution because we have imposed the condition that the third and 
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four moments be satisfied. We would like to say that we have found the Green's function for a strongly correlated 
system, keeping the full (k, oS) - dependence. With our rather simple approach we have been able to have a Green's 
function with cither two poles or three poles. Our approch is somewhat different from the one of the limit of infinite 
dimensions p0[ where only the dynamical properties are taken into account leaving aside the study of the long range 
behavior. The moment approach is a reliable tool to study strongly correlated electronic systems, in particular, the 
Hubbard model. A recent calculation by Nolting, Jaya and Rex [Ell has applied it to the periodic Anderson model, 
where the relevant quantity of study is the self-energy. 
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Figures 



Fig. 1. The self-consistent energy spectra, u>i(k), i = 1,2,3, as function of momentum along the diagonal of the 
Brillouin zone for U/t = 7.0, p = 0.5, T/t = 0.5. 

Fig. 2. ai(k) vs k, i = 1, 2, 3, with the same parameter of Fig. 1. 
Fig. 3. Qi(k) vs k, i = 1, 2 for the same parameters of Fig. 1. 
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